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We revisit the light pseudoscalar Ai in the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) with partial universality at some high unification scale in order to delineate the parameter 
space regions consistent with up-to-date theoretical and experimental constraints and examine to 
what extent this state can be probed by the Large Hadron Collider (LHC) during Run 2. We find 
that it can be accessible through a variety of signatures proceeding via Ai and/or bb, the 

former assuming hadronic decays and the latter two 5-tags within a fat jet or two separate slim ones. 
Herein, the light pseudoscalar state is produced from a heavy Higgs boson decay in either pairs or 
singly in association with a Z boson (in turn decaying into electrons/muons). 


I. INTRODUCTION 

In this report we analyse in detail some of the processes that yield sizeable event rates and could potentially 
lead to the detection of a light Ai at the LHC with ^/s = 14 TeV within the NMSSM [ 1 ] . We perform parameter 
scans of this scenario with partial universality at the Grand Unification Theory (GUT) scale to find regions 
where a light, ^150 GeV, Ai can be obtained. In these scans we require the mass of the SM-like Higgs state 
discovered at the LHG, henceforth denoted by iLsM, to lie around 125 GeV and its signal rates in the 77 and 
ZZ channels to be consistent with the SM expectations. We study in detail the two possibilities, Hsu = Hi 
and Hsm = H 2 , as two separate cases. (Recall that the neutral Higgs spectrum of the NMSSM includes three 
GP-even states. Hi 2 ^ 3 , and two GP-odd ones, Ai^ 2 ^ wherein an increasing numerical label represents an heavier 
state.) Moreover, we assume the Ai to be produced via the decay of a heavy scalar Higgs boson of the model, 
the latter induced by gg fusion. (As finally established in [ 2 ], although some scope was demonstrated for single 
Ai production in association with a bb prior to LHG Higgs data [3], this channel no longer carries any promise. 
Also, note that the scope of Vector Boson Fusion (VBF) and Higgs-strahlung is also currently being re-assessed 
in the light of the same experimental results [4].) In particular, we include the two intermediate channels AiAi 
and AiZ while the decaying heavier Higgs boson can be any of the three neutral scalars. Ai’s thus produced 
decay into either bb or (fully hadronic) r+r“ pairs. The former decay is always the dominant one as the ratio 
of the branching ratios (BRs) for these modes is given approximately by the ratio of the b and the r masses 
squared, but the latter decay can be equally important due to a relatively smaller background. In case 

of the AiZ decay channel, we only consider the leptonic (e+e“ and /i^/i“) decays of the Z boson. 

To study the prospects for the discovery of an Ai at the LHG in all these production and decay channels, 
we employ hadron level Monte Garlo (MG) simulations. We perform a detailed signal-to-background analysis 
for each process of interest, employing jet substructure methods for detecting the b quarks originating from an 
Ai decay, assuming two 5-tags in either two single 5-jets or one fat 5-jet. In particular, notice that, in case of a 
decaying SM-like Higgs state, the mass measurement of ^ 125 GeV serves as an important kinematical handle 
for all signatures. Removing this condition reduces the sensitivity by a factor of 2 to 3. To recap, the AiAi 
pair thus produced decays into the 5555 (45), bbT^T~ (252r) and r+r“r+r“ (4r) final state combinations while 
in the case of AiZ production we will be looking at bM~^i~ {2b2i) and ( 2 r 2 ^), wherein i = e,/a. 


II. MODEL SETUP AND SCANS 

In order to remedy the proliferation of parameters typical of realistic models of Supersymmetry (SUSY), one 
usually invokes some kind of unification of these at high energy scales, typical of GUTs. However, as noted 
in [ 5 ], the fully constrained NMSSM (where all scalar and fermion masses as well as dimensionful couplings 
are unified, respectively, into three separate parameters) struggles to achieve the correct mass for the assumed 
SM-like Higgs boson, particularly in the presence of the latest theoretical and experimental constraints. In order 
to bypass this, the strict unification conditions mentioned above need to be relaxed. In a partially unconstrained 
version of the NMSSM the soft masses of the Higgs fields, 5 taken as independent (from 
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TABLE I: The CNMSSM-NUHM input parameters and their scanned ranges. 


Parameter 

mo (GeV) 

mi /2 (GeV) 

Ao (GeV) 

Meff (GeV) 

Range 

200 - 2000 

100 - 1000 

-3000 - 0 

100 - 200 

tan 13 

A 

K, 

(GeV) 

(GeV) 

1 - 6 

0.4 - 0.7 

0.01 - 0.7 

-500 - 500 

-500 - 500 


mo) parameters at the GUT scale. Through the minimisation conditions of the Higgs potential these three soft 
masses can then be traded at the Electro-Weak (EW) scale for the parameters n, /ieff and tan^^. Similarly, the 
soft trilinear coupling parameters Ax and though still input at the GUT scale, are also taken as independent 
(from Aq). The model is thus defined in terms of the following nine continuous input parameters: 

mo, mi/ 2 , Ao, tan/3, A, n, /ieff, Aa, A^, 

where tan/3 = Vu/vd^ with Vu being the Vacuum Expectation Value (VEV) of the mtype Higgs doublet and 
Vd that of the d-type one. This version of the model serves as a good approximation of the most general EW- 
scale NMSSM as far as the Higgs sector dynamics is concerned. We, therefore, adopt it here to analyse the 
phenomenology of the light pseudoscalar and we refer to it as the GNMSSM-NUHM, where NUHM stands for 
Non-Universal (soft) Higgs Masses. 

We scanned the GNMSSM-NUHM parameter space given in Tab.j^to search for points giving m^^ < 150 GeV. 
We used the publicly available package NMSSMTools-v4.2.1 [ 6 ] for computation of the mass, coupling and BR 
spectrum of the Higgs bosons for each model point. In our scans we imposed the following constraints from 
6 -physics, based on [7], and dark matter relic density measurements, from [8], as 

• BR(Bs ^ /^^/i~) = (3.2 (± 10 % theoetical error) ± 1.35) x 10 “^, 

• BR (Bu ^ ru) = (1.66 ± 0.66 ± 0.38) x 10“^, 

• BR (B ^ Xs 7 ) = (3.43 ± 0.22 ± 0.21) x 10“^, 

• < 0.131 (0.119 ± 10% theoretical error). 

Exclusion limits from the LEP and LHG Higgs boson searches were also tested against using the HiggsBounds- 
v4.1.3 [9] package. Finally, from NMSSMTools we obtained the signal rates of the Hsu state, defined for a 
given decay channel X as 


^ cj{gg ^ X BR{H, ^ V) 

^ ^{99 hsu) X BR(6,sm ^ ’ 

where hsu is the true SM Higgs boson, which we required (for X = 77 , ZZ) to lie within the measured ±lcr 
ranges of the corresponding experimental quantities 

jdjj = 1.13 ± 0.24 , juzz — 1-0 ± 0.29 and (2) 

flryry = l.57^Q 28 5 O'Z Z = ’ (^) 

provided by the GMS [10] and ATLAS [TT] collaborations, respectively. The red and blue points in the 
forthcoming figures are the ones for which the calculated Rx lies within the range of fix measured by GMS and 
ATLAS, respectively, while the green points are the ‘unfiltered’ ones for which neither of these two constraints 
are satisfied. 


III. SIGNAL-TO-BACKGROUND ANALYSIS 

Following the scans, we carried out a dedicated signal-to-background analysis based on MG event generation 
for pp collisions at 14TeV and variable integrated luminosity, for each process of interest. Using the program 
SuSHi-vl. 1.1 [ 12 ], we first calculated the gg fusion production cross section of a SM Higgs boson with the 
same mass as as that of our Hi. This cross section was then rescaled using the ggHi reduced coupling in the 












Toyama International Workshop on Higgs as a Probe of New Physies 2015, 11-15, February, 2015 


3 




[GeV] [GeV] 


FIG. 1: Expected experimental sensitivities as functions of ttiai in various possible final state combinations for (a) 
gg —^ ^ and (b) gg —^ Hi —y A\Z. 


NMSSM and multiplied by the relevant BRs of the Hf^, all of which are obtained from NMSSMTools. The 
backgrounds, which include pp 46, pp 262r, pp ^ At, pp ^ Z2b and pp Z2 t, were computed with 
MadGraph5_aMC@NLO [13]. Both signal and background for each process were hadronised and fragmented 
using Pythia 8.180 m interfaced with Fast Jet-v3.0.6 m for jet clustering. The parton-level acceptance cuts 
used in the event generation in MadGraph are: (i) \r]\ < 2.5 and pt > 15 GeV for all final state objects; (ii) 
AR = ^(Ar/)^ + (^0)^ > 0-2 for all 6-quark pairs; (iii) AR >0.4 for all other pairs of final state objects (where 
Pt^ rj, (j) are the transverse momentum, pseudorapidity and azimuthal angle, respectively). 

As intimated, our use of the jet substructure method m implies that we have three possible signatures 
for a decaying Ai: one fat jet, two single 6-jets and two r-jets. The fat jet analysis, which assumes boosted 
6-quarks, allows one to obtain much higher sensitivities, particularly for large masses of the decaying Higgs 
bosons. Notice that a key ingredient of this selection is the retention of two 6-tags in both cases of a single fat 
6-jet and two slim 6-jets. Failing this, i.e., if only one 6-jet were to be tagged instead, the list of backgrounds 
would dramatically increase, in the form of QCD processes also including light-quark and gluon jets. 

In general, we see in Fig. [^b) that the fat jet analysis can be very effective when the Ai is much lighter 
than the Hi {i = 1,2,3), but gets worse as increases and, in fact, soon becomes relatively useless (the 
corresponding curves are thus cut off at the mass above which the analysis becomes ineffective). This is due 
to the fact that the fat jet analysis assumes boosted 6-quark pairs. One can also see (especially in the curve 
with tuh' = 350 GeV; hereafter H" refers to any of the three CP-even Higgs state directly produced in gg 
fusion while H' refers to the two states other than the FJsm in a given case) that, if the Ai mass becomes 
too small compared to the Hi mass, the sensitivity diminishes due to the 6-jets becoming too collinear to be 
separable even with jet substructure methods. In the upper end, the cut-offs (for sensitivity curves other than 
those relying on the fat jet analysis) are determined by the kinematical upper limit for the given channel, i.e., 
ruAi ~ 62.5 GeV for FJsm ^ and tua^ ~ 35 GeV for FJsm ^ AiZ. 

We finally calculated the expected cross sections for the signal processes which yield S/\/B > 5 for three 
benchmark accumulated luminosities at the LHC, C = 30/fb,300/fb and 3000/fb, in various final state combi¬ 
nations, as functions of Notice that, in order to keep the figures readable, in the following section we will 

only show the curves corresponding to the analyses with the highest sensitivities for a given channel. 


IV. RESULTS 

In the NMSSM, Hi and H 2 can both have masses around 125 GeV and SM-like properties, thereby alterna¬ 
tively playing the role of iLsM- A SM-like Hi with mass around 125 GeV can be obtained over wide regions of the 
CNMSSM-NUHM parameter space defined above. However, the additional requirement of tra^ ^ 150 GeV sig¬ 
nificantly changes this picture. In Fig.j^a) we show the distribution of the mass of Hi against that of Ai for the 
points obtained in our scans assuming Hsu = Hi. We allow a rather wide range of tuhsm^ f22 GeV — 129 GeV, 
in order to take into account the experimental as well possibly large theoretical uncertainties in its model pre¬ 
diction. The heat map in the figure corresponds to the parameter tan/3. One can see a particularly dense 
population of points for tan ^ 1 — 6 in the figure, with the mass of Hi reaching comparatively larger values 

than elsewhere. However, itiai for such points almost never falls below ~ 60 GeV. In Fig. [^b) we show tuhi 
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FIG. 2: Case with i/sM = Hi: (a) Mass of Hi vs. that of Ai, with the heat map showing the distribution of tan/3; (b) 
mHi as a function of the parameter n, with the heat map showing the distribution of the coupling A. 



m^JGeV] 





FIG. 3: Mass of H 2 vs. that of Ai for the case with Hqm = H 2 . The heat map shows (a) the distribution of tan/3 and 
(b) the distribution of A. 


as a function of the coupling k, with the heat map corresponding to the coupling A. Again there is a clear strip 
of points with A ^ 0.6 (and n ^ 0.15 — 0.5) for which ruHi can be as high as 129 GeV. The rest of the points, 
corresponding to smaller A and larger tan/3, can barely yield tuh^ in excess of 126GeV. The reason for the 
behaviour of tuh^ observed in these figures is well explained in Ref. [2] and is essentially attributable to the 
expression of in terms of the CNMSSM-NUHM input parameters. 

Figs, [^a) and (b) show that H 2 with a mass lying in the entire allowed range can be obtained much more 
easily without always requiring very low tan /3 or very large A. Moreover, the corresponding parameter space 
points can also yield fairly small Ai (with sizeable BR(i72 ^ AiAi/Z)), without the H 2 deviating too much 
from the LHC Higgs boson signal rate measurements. We will, therefore, concentrate in the remainder of this 
report only on the H 2 solution for iLsM (i-c., i^sM = H 2 ). 


A. Production via Hqm AiAijZ 

In Fig. [^a) we show the prospects for the H 2 AiAi channel when H 2 is SM-like. We see that, compared 
with the iLsM = Hi case, a much larger part of the parameter space can be probed at the LHC, even at as low 
as 30/fb of integrated luminosity. The reason is clearly that in this case the points with tuai < belong 

to the parameter space regions where BR(iL 2 ^ AiAi) is indeed sufficiently enhanced without causing the Hsu 
for these points to depart from a SM-like behaviour. This is also the reason why a large fraction of the points 
with large event rates is consistent also with the CMS and ATLAS measurements of Fig. j^b) we 

see instead that the prospects in the H 2 AiZ channel are poor. 
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FIG. 4: Total cross sections for i/sM = H 2 for (a) gg Hqm AiAi and (b) gg Hqm AiZ. 
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FIG. 5: Total cross sections for i/sM = H 2 for (a) gg ^ Hi ^ AiAi and (b) gg ^ Hi ^ AiZ. 


B. Production via H' AiAijZ 


The prospects for the discovery of a light pseudoscalar in the Hi —> AiAi and Hi AiZ decay channels, 
for a singlet-like Hi, are illustrated in Figs, ^a) and (b), respectively. One sees in Fig.j^a) that almost all the 
points complying with the current CMS and/or ATLAS constraints on Rx are potentially discoverable, even at 


jC = 30/fb. Thus a large part of the scanned NMSSM parameter space can be probed via this decay channel. 
In particular, since such light pseudoscalars cannot easily be obtained for the case with Hsu = Hi, it should 
essentially be possible to exclude or confirm < 60 GeV in the NMSSM at the LHC via this channel. Note 
also that such an exclusion will not cover the narrow regions of parameter space where trai > Finally, 

In Fig. ^b), we see that the prospects for the discovery of Ai via the Hi AiZ channel are non-existent. 

For the decay chain starting from H^, the situation is illustrated in Fig.j^a), where we see that the H^ AiAi 
channel is inaccessible also due to the fact that, for such high masses of H^ 400 GeV), the production cross 
section gets diminished. Moreover, other decay channels of Hs dominate. Conversely, the Hs AiZ channel, 
shown in Fig. |^b), shows much more promise. This has to do with the increased sensitivity in the fat jet 
analysis when the involved masses are high as well as the relatively large H^AiZ coupling, which is actually 
somewhat larger here than in the Hsu = Hi case, due to a correspondingly larger doublet component of Ai. 
We therefore emphasise again that this channel will be an extremely important probe for an NMSSM Ai with 
mass greater than ^ 60 GeV. 


V. SUMMARY 

We have found that the decays of the NMSSM CP-even scalars, including in particular the SM-like Higgs 
boson, whether Hi or H 2 , carry the potential to reveal an Ai with mass < 60 GeV for an integrated luminosity 
at the LHC as low as 30/fb. This is particularly true when the SM-like Higgs state is the H 2 . Most notably 
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FIG. 6: Total cross sections for i/sM = H 2 for (a) gg H 3 AiAi and (b) gg H 3 AiZ. 


though, when the Ai is heavier than ^ 60 GeV, while its pair production via decays of the two lightest CP-even 
Higgs bosons becomes inaccessible, the gg AiZ channel takes over as the most promising one. This 

(hitherto neglected) mode is, therefore, of great importance and warrants dedicated probes in future analyses 
at the LHC. 

Such results are based upon parton-level MC event generation supplemented by parton shower and hadroni- 
sation. Hence, final validation of our results can be done through a proper detector simulation. A key aspect 
of this would be the implementation of double 6 -tagging throughout. 
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